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Abstract 

A relevant relation between the dwell time and the density of states for a three 
dimensional system of arbitrary shape with an arbitrary number of incoming 
channel is derived. This result extends the one obtained by Gasparian et al. 
for the special case of a layered one dimensional symmetrical system. We 
believe that such a strong relation between the most widely accepted time 
related to the dynamics of a particle and the density of states in the barrier 
region, one of the most relevant properties of a system in equilibrium, is rich 
of physical significance. 

PACS numbers: 03.65.-w, 71.20.-b, 73.20.At, 73.40.Gk 
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I. INTRODUCTION 



In the controversial field of tunneling times, the dwell time is widely accepted as the 
average time spent by a particle in a given region of space. This time was postulated 
by Biittiker in the context of the tunneling time problemS; recently, a rigorous derivation 
of the dwell time was obtained within Feynman'fl and Bohm'si formulations of quantum 
mechanics. We should also mention that the agreement on the physical significance of the 
dwell time is not unanimousiH. 

The connection between the dwell time and the density of states in the barrier region 
was shown by Gasparian et al.iS, for the case of a one dimensional symmetrical multilayered 
system. In this paper we will extend this result to an arbitrary three dimensional region, 
using a very concise derivation. 

We believe that establishing and clarifying the connection between these two quantities 
is also relevant to the recent work by Wang et al.i, whose object is the investigation of 
the statistics of quasi-bound states in classically chaotic systems. These authors assume a 
relation between the energies of quasi-bound states and the extrema of the dwell time in the 
system under consideration, which can be addressed using the results we present here. 



II. DWELL TIME AND DENSITY OF STATES 

Let us consider a region Q in three dimensional space, connected with the outside world 
by channels. Let a given channel be characterized not only by its spatial location, but 
also by its particular propagation mode. 

We choose as an orthonormal basis for this system the set of stationary state eigenfunc- 
tions corresponding to each of these incoming channels, i.e. {\(pn{E))} , where n = 1 . . . N. 
Our basis is continuous and degenerate, so that an appropriate normalization is 

{ME'MmiE")) = SnmS{E' - E")' (l) 

of course, we also have 

^ r 

Y: dE\UE)){ME)\ = l. (2) 

n=l 

We now derive an expression for the dwell time associated with a general wave function 
in this set. Let us consider a wave packet \tpn{t)) incoming from the nth channel, such that 
the probability of finding the particle in Q vanishes for time approaching ±00. can 
be written in terms of our basis as 

\Mt)) = J an{E)e-'^'/^\ME))dE, (3) 

where / \an{E)\'^dE = 1 so that lipnit)) is normalized to unity. 

The mean dwell time in Q associated with the wave packet i 



r 



00 Jo, 
00 



dt{i;4t)\Pn\Mt)), (4) 
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where Pq in the projection operator onto the region n. Substitution of (D in §) yields 



ri^^ = 2Tin J \an{E)\^{ME)\Pn\ME))dE. (5) 

Therefore the dwell time t^\e) for the stationary state \(j)n{E)) can be defined as 

tP{E) = 27rn{cPn{E)\Pn\ME)), (6) 

which corresponds to the limit of @ as tends to a delta function. 

Equation (H) is equal to the well-known expression postulated by BiittikeiS, provided the 
different normalization for the wave function used here is taken into account. In fact, for 
the wave functions of our basis the incoming probability current is just (27r^)~^ (the total 
probability is given in units of inverse energy, according to normalization (|I|), therefore the 
probability current is in units of inverse action). In the Appendix two interesting formu- 
las relating the dwell time to perturbative potential approachesiH and to the Green's 
functions for our system are shown. 

The local density of states p(r, E) is given by0 

N 

p{r,E) = Y: / {r\UE')){ME')\r)5{E-E')dE', (7) 

n=l 

which, in our case, becomes 

TV 

p{v,E) = J2{r\UE)){UE)\r). (8) 



n=l 



The density of states pn{E) for the region fl is just the integral of p(r, E) over Q, therefore 
we obtain 



TV 

Pn{E) = Y: / {r\ME)){ME)\r)dr 



N 



j:{ME)\Pn\ME)). (9) 



n=l 



From d^) and (^, one straightforwardly obtains 

l'^(E) = —Y.r'S\E), (10) 

n=l 

i.e., the density of states in ^2 is proportional to the sum of the dwell times in Vl for all 
incoming channels. 



III. COMMENTS 



The result of Gasparian et al.i can be easily obtained as a particular case of (0). In 
fact, for a one dimensional region the number of channels reduces to two (for left and right 
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incoming waves) and for a symmetric potential we have Tjj\E) = t^\e) = td{E), so that 
we can write 

^[Idim-sym] = (^^) 

which corresponds exactly to Equation (5) of Ref. ^. 

If the region Q is connected to the outside world by tunneling barriers, so that we 
have quasi-bound states in Q, the density of states is strongly peaked for energy values 
corresponding to these quasi bound states (in the limit of a closed region, the density of 
states has to be a set of delta functions). Searching for the peaks in the density of states is 
therefore a possible way to find the statistics of quasi bound states. Moreover, formula (|TUp 



tells us that a peak in density of states corresponds to a peak in the dwell time for one of 
the incoming channels, so that, in order to evaluate the statistics of quasi bound state Jl, it 
is correct to search for the maxima of the dwell times of all the possible incoming channels. 
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APPENDIX A: TWO FORMULAS FOR THE DWELL TIME OF STATIONARY 

STATES 



In this Appendix we want to demonstrate two formulas which relate the expression (^ 
for the dwell time to perturbative potential approacheJlEi and to the Green's functions for 
the system under consideration. 

Let us apply a perturbative potential V to the region Q, so that VPq is the perturbation 
operator to be added to the unperturbed Hamiltonian H. The new orthonormal basis for 
this system is {{(j)^ y{E))} , where each function is a solution of the equation 



{E±e-H-VPnMtvm = 0, 



(Al) 



for e — > and the nth incoming channel. 

We want to show that the dwell time expression given by (P) is just the diagonal matrix 

0, i.e., 



element of 2ih-^ evaluated for V 



tP{E) = [2th 



-] 



v=o 

{ct>iy{E)\2thi,\c^Xy{E)) 



{<P^yiEMyiE)) 



(A2) 



v=o 
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Let us point out that both the numerator and the denominator diverge (because the basis 
is Dirac-normahzed) . Nevertheless their ratio is finite, as we will show. 
We can write0 



\<P%ym = ± VG^{E)Pa\4>lvm. (A3) 

where G^{E) is a solution of the equation 

{E±e- H)G^{E) = 1, (A4) 
for e ^ 0. Substituting in yields 

(0„(E)|0„(E)) 

Since in the righthand side of this equation is a diagonal matrix element, it has to be real, 
therefore only the imaginary part of G'^lE) matters. But we know0 that 

N 

lm{G^{E)} = Tvr E (A6) 

m=l 



which, substituted in (|A5|) allows us to write 

[E) = 2Tin 2^ 

= 2ixh{<Pn{E)\Pn\<pn{E)) (A7) 
where the last step derives from the orthogonality of the basis. Therefore Equation (^) is 



shown to be equivalent to (|A2|) and (|A5|). 
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